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Abstract 

In this paper we present a systematic view of Quan- 
tum Cellular Automata (QCA), a mathematical for- 
malism of quantum computation. First we give a 
general mathematical framework with which to study 
QCA models. Then we present four different QCA 
models, and compare them. One model we discuss is 
the traditional QCA, similar to those introduced by 
Shumacher and Werner, Watrous, and Van Dam. We 
discuss also Margolus QCA, also discussed by Schu- 
macher and Werner. We introduce two new models, 
Coloured QCA, and Continuous-Time QCA. We also 
compare our models with the established models. We 
give proofs of computational equivalence for several of 
these models. We show the strengths of each model, 
and provide examples of how our models can be use- 
ful to come up with algorithms, and implement them 
in real- world physical devices. 
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Figure 1: Wolfram's "Rule 30": Each cell has a bi- 
nary state, and it is updated according to the value 
of itself, and its two nearest neighbours, one to each 
side. There are eight possible states for the neigh- 
bourhood of any particular cell. The diagram above 
gives the possible neighbourhood configurations, with 
the corresponding updated state for the centre cell 
below. Reading the second row as binary digits, one 
obtains Wolfram's name for this rule: 30. 



previous models, whether given the name of QCA or 
some other, are either instances of, or equivalent to, 
our model of QCA. Along the way, we will survey 
some of the major work in the field of QCA. 



1 Introduction 



2 Classical Cellular automata 



Quantum cellular automata (QCA) research has seen 
significant growth in the recent years. This model 
of computation has been appearing in the literature, 
sometimes with different names, and in several dif- 
ferent guises. Quantum lattice gases, pulse-driven 
quantum computers, and translation-invariant quan- 
tum operators, are all instances of QCA, and yet little 
has been discussed of the relationship between these. 

The purpose of this paper is to give one unifying 
view of all QCA type phenomena, in the form of a 
model of computation. We intend to show how all 



We start our discussion with a short review of classi- 
cal cellular automata. 

A Cellular Automaton (CA) consists of a lattice 
structure, where each cell is in one of a finite number 
of predetermined cell states. At each discrete time- 
step, every cell is updated, in parallel, according to a 
local, spatially uniform rule. This gives us a model of 
computation which is different from the usual circuit 
model or the Turing machine model. 

An example of a CA rule is given in figure ^ I n 
brief, a CA is a state lattice with an update rule that 
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Figure 2: Rule 30 in motion: Here we can see several 
steps of the CA Rule 30 as presented in Figure H 
This is a one-dimensional lattice, with the downward 
vertical dimension representing successive time steps. 

updates the value of each lattice cell at each time 
step. The rule is local, in that the updated state for 
an individual cell depends only on the current states 
of the cell itself and of its neighbour cells. The set of 
cells on which the updated value of cell x depends is 
called the neighbourhood of x. The size and form of 
the neighbourhood may vary from CA to CA. 

We can see the time evolution of the rule in Figure 

□ 

In order to quantize the CA model, it is desirable 
to first make the model reversible. On the other 
hand, reversible computing predates quantum com- 
puting by several decades. One way to ensure that 
a CA is reversible is to use a Margolus partitioning 
scheme, in which the lattice is divided into tiles. A 
reversible operation is performed on each tile. On 
consecutive time-steps the tiling is staggered to allow 
the possibility of data propagation (see Figure |3J). 

2.1 Quantum Cellular Automata 

We wish to to formulate a model for quantum cellular 
automata. Several properties are desirable: 

• The model should be a generalization of the clas- 
sical CA and subsume the latter as a special case. 

• The model should allow for quantum computa- 
tion; it should be equivalent to the quantum cir- 
cuit model. 
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Figure 3: Margolus Partitioning CA Scheme 

• The model should ideally be reversible. 

• The model should be a natural abstraction for 
quantum computation on particular realistic de- 
vices. 

With this in mind, we first present a generic, basic 
model which will be our starting point for all our 
QCA studies. 

2.1.1 The Basic Model 

In our quantum cellular automaton model, we con- 
sider qubits arranged in an integer lattice of dimen- 
sion d, Z d . 

Instead of trying to define an state on an in- 
finite lattice, we define a state p as a family of 
states on finite subsets of the lattice. For each fi- 
nite A C Z d , we take pa to be a state in the Hubert 
space Tt® A over the subset A. The family of states 
p = {pa} Acz d ,A finite must also satisfy a consistency 
condition. For every finite subset A C Z d , given 
B C. A, we must have 

Pb = ^ A \b{pa)- 

In other words, partial traces must be consistent. 

For a QCA, we must also have a global evolution 
operator. However, for any particular qubit in the 
lattice, interactions in one time-step must be local- 
ized to within a given neighbourhood of the qubit. To 
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this end, we will first define a neighbourhood scheme 
as a finite set N of elements of Z d which includes the 
zero vector 0, and then for each given lattice posi- 
tion x, we define the the neighbourhood of x to be 
Af(x) = x + N. For a finite subset A C Z d , we will 
define JV(A) to be the union of the neighbourhoods 
of the elements of A. Naturally, neighbourhood sets 
are translation independent, that is, given a lattice 
translation r, we have J\f{r{x)) = t{M{x)). 

Now, given a finite subset A C Z d , the result of 
a global evolution operator can be determined from 
its action on the subset Af{A) . Accordingly, we define 
our global evolution operator U as a family of unitary 
operators Ujj-^a) acting on neighbourhoods for every 
finite subset A C 7L d . The action of U on a given 
state p is given by 

{Up) A = ^N{A)\A {UN(A)PN(A)Ul /{A ^j . 

Of course, we also require U to satisfy translation 
independence, so that for any subset A C Z d and 
lattice translation r, we have t//v"( r (A)) = t {UaT(A))- 

2.1.2 Deriving Global and Local Operators 

We consider the operator U, consisting of a unitary 
operator Uj^-(a) for each finite subset A C Z d , as the 
global evolution operator for a given cellular automa- 
ton. From this, it would be desirable to find a local 
operator, which could reciprocally be used to derive 
the global operator. 

The unitary operator Uj^^ , corresponding to the 
subset {x} consisting of a single element, makes a 
natural choice for a local operator. However, deriving 
the global operator from this requires some work, as 
we require a family of unitary operators which satisfy 
the consistency condition for each finite subset of Z d . 

To this end, for the qubit at lattice position x, we 
define S° to be the set of pure states \tp) over the 
lattice subset J\f{x) which satisfy 

and similarly, we define Si to be the set of pure states 
which map to |1)(1|. In general, the set S 3 X is an affinc 
space for which a basis may be computed for any 



given Ujyr x y We may also consider the affine spaces 
Si as sets of states over a larger finite subset of Z d 
containing x. 

Now, given a finite subset A C Z d , with A = 
{x\ ,x n }, the set of lattice states which maps to a 
particular state \ j) — \ ji, ■ ■ ■ ,j n ) on the lattice subset 
A will be 

n 

s °a = n s *k > 
fc=i 

where the sets Si k are considered as sets of states 
over the lattice subset M{A) . 

By selecting a basis for the set of states over the lat- 
tice points in J\f{A)\A for each j, it is possible to find 
a unitary operator LV(A) that makes the appropriate 
map, so any vector in S A is properly mapped to 
In addition, note that in order for Ua to satisfy the 
consistency condition, vectors outside S A must not 
be mapped to \j). Since every vector in S A must be 
mapped to something, it must be \ j). Thus every op- 
erator U that satisfies the consistency condition can 
be constructed in this way. Also, note that while we 
have many choices for any particular operator Ua, 
every choice yields the same result after the partial 
trace. Next, we will show these operators are always 
consistent. 

In order to be a proper global transition operator, 
U must preserve the consistency condition. That is, 
given a state p, and given that a — Up, we must have 

oa = ^b\a{^b) 

for every finite subset B C 7L d and for any A C B. In 
other words, we must have 

tT AT(A)\A (Uj^^A)PAr(A)Ul r ^ A ^j = 

(W( B)PN-(B)Ul f{B ^j , 

given that trx(B)\N{A)PN{B) = PaT(A) 

Since both Uj^-^a) , C^V( b) an d the partial trace are 
linear functions, it suffices to show that this relation 
holds for operators of the form p//m) = 1^X01 > where 
\ip) € S 3 B and |0) E S B for computational basis states 
\j) and \k) over B. But by definition, we also have 
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Figure 4: Margolus QCA: Similar to the classical re- 
versible CA, the lattice is partitioned using two dif- 
ferent tiling sets. Each tile is then acted on with a 
unitary operator. 

\ip) € S° A and |</>) £ S\, and so both sides of the 
above identity equal trm^|j)(fc|. 

Therefore, an operator U is valid if and only if it is 
the same as the one constructed by the local operator, 
up to unitary operators acting on Af(A) \ A. 

3 QCA Models 

Now that we have shown our general QCA model, 
we can introduce particular varieties of QCA types 
as instances of our model. The first one we introduce 
is the Margolus QCA, or MQCA. 

3.1 Margolus QCA 

The MQCA is similar to its classical counterpart, in 
which we define a tiling of the lattice, and apply the 
same unitary operation to each tile. The tiling is then 
changed in such a way that each tile of one tiling 
overlaps with at least two tiles of the original one. 
The two tilings need not have equally shaped tiles, 
but must have the same period with respect to the 
lattice (see Figure 0}. 
Formally, we state: 

Definition 1. A Margolus partitioning is a pair of 
array partitions k and r, dividing the lattice into col- 



lections of identical, finite, disjoint, and uniformly 
arranged blocks, such that each block of one partition 
overlaps with at least two from the other. A Margo- 
lus Quantum Cellular Automata (MQCA) is defined 
by lattice S a Margolus partitioning of S and a pair 
U, V of unitary transformations, each one acting on 
the blocks of one of the partitionings. 

It is easy to see that the MQCA follows the con- 
sistency conditions outlined above. An interesting 
example of an MQCA is the multi-particle quantum 
walk. 

Example 1. Let S be C z , the integer line. The par- 
titionings p and s both divide the lattice into blocks of 
two qubits, each block overlapping one qubit from one 
block in the other partitioning, as in Figure^ Let 

/ 1 \ 



\ 1 / 

For this example, in the computational basis, we 
understand a 1 to represent the presence of a particle, 
and zero to represent the abscence. The overall effect 
is that any particles in the lattice slowly get diffused 
over time. There is no interaction among particles. 

Both van Dam |TQj and Watrous consider parti- 
tioned QCA or PQCA. A PQCA is shown in detail 
if figure [3J This scheme is simply a particular kind 
of Margolus partitioning; one of the Margolus parti- 
tioning operators is simply a swap. 

It should be pointed out, though, that the van Dam 
and Watrous model differs greatly from the one pre- 
sented here. In their model, the quantum state of 
QCA is in a superposition of classical configurations 
of the lattice. This leads to a very different defini- 
tion of the evolution operator as well, which runs into 
problems as outlined in [3]. Type-II quantum [rollT7| . 
computers can be regarded as particular instances of 
the PQCA scheme, except that the exchange proce- 
dure is done non-coherently. 

The example we chose also has a historical signif- 
icance. Quantum walks on a lattice have been ex- 
tensively studied before (see for example 0] for an 
overview)), and in particular a QCA quantum walk 
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3.2 Colored QCA 
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Figure 5: Partitioned QCA: In a PQCA each cell is 
partitioned into 3 'subcells': left, centre, and right. 
A unitary operator is applied to each cell at odd time 
steps, and at even time steps the cells exchange sub- 
cells, it's right subcell is exchanged with the left sub- 
cell of the right neighbour, and the left subcell is ex- 
changed with the right subcell of the left neighbour. 



has been analyzed in [o]. It is of note that one can 
take the continuous limit of the walk we outlined 
above, i.e. the limit when the space between lat- 
tice cells and the duration of the time step both go 
to zero. Doing so gives the Schrodinger equation for 
a particle in free motion This works for lattices 
of any dimension. By taking the limit in a slightly 
different way one can obtain the Dirac equation; al- 
though this works only for one dimensional systems 
at the moment [H]- 

Also, in 3 Schumacher and Werner show that 
MQCA are universal QCA, in the sense that all QCA 
(as defined in their general scheme) can be reduced 
to MQCA. 

More recently, Raussendorf shows how to build a 
universal quantum computer using a Margolus style 
quantum cellular automata [7]. For these reasons, 
implementing a MQCA on a physical device would 
be of great usefulness. We explore this in our next 
section. 



In 1993 S. Lloyd introduced an example of a (time- 
inhomogenous) CQCAs, which we will call Spin- 
chain QCA. He introduced this model as a proposal 
for a feasible quantum computer. Consider a one- 
dimensional chain of spin 1/2 systems, such as a poly- 
mer, with three different species, i.e ABC ABC . . . 

Suppose that the nearest-neighbour interactions 
are given by (arbitrary) Hamitlonians Hab, Hbc-, 
He a- In the computational basis, the diagonal terms 
of this Hamiltonian shift the energy levels of each cell 
as a function of the energy levels of its neighbours |S] . 

The resonant frequency lo a can then take the value 
u>oq, ljq!, u)i , or u)^, depending on whether its C and 
B neighbours are in the states and 0, and 1, 1 and 
0, or 1 and 1. 

Hence, chain link A with neighbours C and B has 
resonant frequency oja with distinguishable diagonal 
terms Wq , u)^, uj^q, and w^. If these are all differ- 
ent, then transitions on species A spins can be done 
selectively depending on the value k. For instance, by 
applying a n pulse with frequency luq Q all species A 
lattice points whose both neighbours are in the state 
will be flipped. 

It is also possible to apply, using the same tech- 
niques, any arbitrary single-qubit gate to all spins of 
the same species and to apply any two qubit gate on 
all pairs A, B or C, A or B, C. 

We introduce a new general QCA model, which we 
call Coloured QCA (CQCA). This model is a gener- 
alization and abstraction of Lloyd's scheme. 

In a CQCA, each lattice point is assigned a colour 
in a checkerboard fashion. At each time step only 
points of a certain colour are updated with a unitary 
dependant on their neighbour's values. Neighbours 
of the same colour are not distinguishable. 

A graphical representation is shown in Figure H3 
and a formal definition is given below. 

Definition 2. A correct k- colouring for a lattice is a 
periodic mapping C from lattice points to {0, . . . , k — 
1}, such that no two neighbours have the same colour. 
Fix a single qubit observable a. For each lattice cell 
t its I -field is the value J2 r eN(t) c(r)=;('^l <T '' r ' ) l' ! /')- ^ 
field- controlled unitary u is such that it is only applied 
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0-00000 

Figure 6: Colored QCA This is a representation of a 
one-dimensional CQCA, with time going downward. 
The two shades of gray differentiate between the two 
species in this CQCA. The brackets represent the 
controlled unitaries: at each step either the light or 
dark grey sites are acted on with unitary operators 
that depend on their left and right neighbour values. 

to a lattice point if its I- field, for I S {0, . . . , k — 1} 
has a particular set of values. 

A Coloured QCA or CQCA is defined by a neigh- 
bouring N , colouring C , and a field- controlled unitary 
ui for each colour I 6 {0, . . . , k — 1}. 

Again, it is not hard to see that CQCA adhere 
to the consistency constraints outlined in our generic 
model. 

Several proposal for implementation of quantum 
computers can be seen as instances of the CQCA 
model. Obviously, Lloyd's pulse-drvien quantum 
computers are one example. Benjamin's proposal [H] 
is an example of 2-CQCA scheme that achieves uni- 
versality. K. G. H. Vollbrecht and J. I. Cirac [TDJ 
present a 1-CQCA (each lattice site being a qudit 
with d = 4), that also achieves universality. 

Example 2. An example of a proper CQCA is the 
1-D guantum walk CQCA. As in Example 1 above, 
let S be C z , the integer line. The colouring scheme 
uses four colours {ci}o<j<3 such that position Xi on 
the lattice has colour i mod 4. 

The evolution operator is as follows. First, c\ is 
given a tt z-rotation conditioned on cq being in 
Then cq is given a tt/2 z-rotation conditioned on c\ 
being in |1). Then we repeat the ir z-rotation on c\ 



conditioned on cq being in We can see that this 
procedure is equivalent to doing a square-root-of-swap 
on all pairs cq-c\. We do the same procedure on 
colours C2-C3, and then on both pairs C\-c<i and C3- 
c Q . 

The above construction gives the exact dynamics as 
the MQCA of Example 1 above, namely, a quantum 
walk on the line. 

From the example above, one might deduce that 
any dynamic given in the MQCA model can also be 
constructed using CQCA. This is in fact true, and is 
proven below. 

Theorem 1. For every MQCA there exists a CQCA, 
and for every CQCA there exists and MQCA that has 
the exact same dynamics. 

Proof. For this proof we rely on the fact that arbi- 
trary single qubit operations, coupled with c-NOT 
gates on adjacent qubits (for a connected topology) 
is universal for quantum computation. 

We proceed constructively. Given an MQCA M, 
with partitions s, t, we can build a CQCA C that has 
the exact same dynamics as follows. 

The lattice of M will be the same lattice for S. 
Without loss of generality suppose that the partition 
s has more lattice sites per block than t. Let k be 
the number of sites in each block of s. We construct 
our CQCA with 2s colours. The colour mapping is 
such that each partition block of M is represented by 
s distinct colours in C (each block having the exact 
same colours as each other block in the partition), 
and two neighbouring block has two similar colours. 

Now, we can simulate any arbitrary unitary acting 
on a block of M using only colour-controlled opera- 
tions in C: since each spin in the block has a dis- 
tinct colour, they are all individually addressable. 
Moreover, we can apply c-NOT gates between any 
two neighbouring spins. Operations on C will be re- 
peated periodically over the lattice, due to the re- 
peating colour scheme, and hence it is important that 
block boundaries coincide with colour periodicity. 

Also we use 2s colours instead s, so it is necessary 
to repeat the operations on C for alternating blocks 
of M. This is necessary to ensure that we can isolate 
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Figure 7: Continuous Time QCA: We have a crys- 
tal with two types of nuclei A and B, one represented 
as light gray spheres, the other dark gray. The lines 
connecting the spheres represent the nearest neigh- 
bour couplings, with different colour lines represent- 
ing different coupling Hamiltonians. Notice that the 
coupling Hamiltonian depends only on the colour of 
the coupled spins. 

each block in order to perform arbitrary operations 
on it. 

□ 

The importance of the above theorem, beyond 
mathematical curiousity, is that it allows us a sim- 
ple way to implement MQCA algorithms on physical 
devices. 

3.3 Continuous-Time QCA 

Another model we present is the Continuous Time 
QCA (CTQCA). 

A CTQCA is similar to a CQCA in that all lat- 
tice points are coloured. However, instead of unitary 
operators applied in discrete timesteps, the system 
evolves continuously according to a Hamiltonian de- 
scribed only by nearest neighbour couplings. 

Definition 3. Let C be any colouring, that is any 
periodic mapping from lattice points to {0, . . . , k— 1}, 
but not necessarily 'correct' in the sense of Definition 
2 above. 

Every pair of neighbouring lattice points i,j has 
a coupling Hamiltonian Hij(t) that depends only on 



the colour of the two points. For a given region of the 
lattice, the Hamiltonian of that region is simply 

The evolution of the CTQCA over a given time 
period t is given by 




A good example of a CTQCA is the diffusion au- 
tomata as defined below. 

Example 3. Again, as in example 1 above Let S be 
C z , the integer line. The colouring scheme requires 
only one colour, and each spin has only two neigh- 
bours, one directly to its left, and to its right. The 
coupling Hamiltonian is 

i 

In solid-state NMR this is called the flip-flop cou- 
pling, since it has the effect of flipping two contiguous 
spins from 1 01) to |10). 

The evolution of this CTQCA is similar to the 
quantum walks of Examples 1 and 2 above. Again 
this is not coincidental. 

An important result is that CTQCA is, again, an 
equivalent computational model to the other QCA 
models presented above. 

Theorem 2. Given, a CQCA there is a CTQCA that 
has the same dynamics. 

Proof. This construction is very simple. Given a a 
CQCA C we can construct a CTQCA T with the 
same lattice and same colour-scheme. We give a time- 
dependant Hamiltonian H that changes at discrete 
time-steps of duration t so that U — e ltH where U is 
the unitary of the CQCA at the given timestep. 

□ 

The converse is slightly more complicated, and can 
only be done in an approximate sense, 
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Theorem 3. Let AT be the time it takes to do one 
time step. Then, a CQCA C can approximate a 
CTQCA T when AT — > 0. 

Proof. Continuous time QCA, as opposed to CQCA 
do not have the restriction that neighbours cannot be 
of the same colour. Assuming that the colouring is 
correct for the given CTQCA T, then we proceed as 
follows. 

Let C have the same colour scheme as T. 

For any pair of colours U and t j let H%j be the cou- 
pling Hamiltonian for neighbours of that colour. Let 
Uij[t] — e~ ltH . Now, Uij[t] acts only on two spins, 
and hence can be described as a series of controlled 
operations, uff ] \t], u\f\t], u\f ] [t], U^f ] \t\... 

where U- j [t] is an operation on the spin of colour i 

controlled by the spin j, and so forth. These U-j S \t] 
become the controlled operators in our CQCA C. By 
letting t <— the evolution of C approximates the 
evolution of the CTQCA T. 

If T has neighbours of the same colour, we simply 
break any colour that has neighbouring spins into two 
or more colours, such that there are no longer neigh- 
bouring spins of the same colour. All the colours 
created this way will have the same coupling Hamil- 
tonian to each other as the original colour had with it- 
self. We can then apply the procedure outlined above. 

□ 



4 Conclusions 

In closing we wish to give a summary of the contribu- 
tions of this paper. First we gave a simple scheme for 
QCA that is general enough to encompass previous 
models of QCA, as well as other phenomenae studied 
under different names. 

We gave several specific models, under this scheme, 
and we showed the relationship among this models. 
We posit that our model is not overly general, that 
is it describes all well formed QCA style phenomana 
and nothing else. This latter statement, though, is 
posited without proof. 
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